Geometric characterization of separability and entanglement in pure 
Gaussian states by single-mode unitary operations 
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We present a geometric approach to the characterization of separability and entanglement in 
pure Gaussian states of an arbitrary number of modes. The analysis is performed adapting to 
continuous variables a formalism based on single subsystem unitary transformations that has been 
recently introduced to characterize separability and entanglement in pure states of qubits and qutrits 
arXiv:0706.1561 . In analogy with the finite-dimensional case, we demonstrate that the 1 x M 
bipartite entanglement of a multimode pure Gaussian state can be quantified by the minimum 
squared Euclidean distance between the state itself and the set of states obtained by transforming 
it via suitable local symplectic (unitary) operations. This minimum distance, corresponding to a, 
uniquely determined, extremal local operation, defines a novel entanglement monotone equivalent to 
the entropy of entanglement, and amenable to direct experimental measurement with linear optical 
schemes. 

PACS numbers: 03.65.Ud, 03.67.Mn 



I. INTRODUCTION 



The concept of entanglement has gradually developed 
from the status of a puzzling interpretational problem, to 
that of a crucial operational resource for quantum infor- 
mation tasks and, even more remarkably, to the status 
of a founding property of quantum theory, whose im- 
plications and applications extend into many diverse ar- 
eas of research ranging from quantum optics and atomic 
and molecular physics to condensed matter physics and 
quantum critical phenomena P, While many open 
questions, even on defining grounds, stand open when it 
comes to address questions like the nature of multipar- 
tite entanglement and the entanglement of mixed states, 
a fairly satisfactory classification and quantification of bi- 
partite entanglement of pure quantum states have been 
established [JHI- This achievement has been partly pos- 
sible because the milestones of quantum information sci- 
ence, like quantum teleportation, quantum cryptography, 
state transfer, broadcasting and telecloning, entangle- 
ment creation and distillation, all rely on the paradigm 
of two distant labs operated by two parties - tradition- 
ally named Alice and Bob - who wish to communicate, 
possibly sharing a pure entangled state [![. By properly 
defining figures of merit associated with such protocols, 
pure-state bipartite entanglement has been understood 
both qualitatively - entangled states are non-separable 
- and quantitatively - the degree of pure-state bipartite 
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entanglement is equal to the entropy of the reduced state 
of each subsystem. In particular, the von Neumann en- 
tropy of entanglement is equal both to the, operationally 
defined, distillable entanglement and entanglement cost 
of pure bipartite quantum states [f| • This equivalence is 
lost in the presence of mixedness, and the phenomenon 
of entanglement conversion irreversibility sets in Q ■ 

There exists an alternative, equally natural way to un- 
derstand and characterize entanglement. From a geomet- 
ric perspective, the degree of entanglement in a state p 
can be quantified as the minimum distance, suitably mea- 
sured, between p and the set of unentangled, separable 
states Q. Again, for pure states, if such distance is mea- 
sured in terms of the relative entropy, the resulting mea- 
sure of entanglement coincides with the von Neumann 
entropy of entanglement Q. This suggests that other 
entanglement monotones, that can be useful either for 
their operational meaning and/or for their computabil- 
ity, might be endowed with an alternative, geometric in- 
terpretation which adds to their understanding and may 
provide alternative tools in their experimental detection. 
A novel approach to the study of this problem has been 
recently introduced for low-dimensional discrete- variable 
systems such as qubits and qutrits [|[. It relies on the 
basic idea that entanglement can be characterized by the 
response of a system to local and unitary perturbations 
that, by definition, cannot change the degree of entangle- 
ment present in the system. Notwithstanding this simple 
fact, oddly enough, degrees of freedoms that are affected 
by local unitaries tend to be systematically neglected in 
the analyses of entanglement properties |{|. The recent 
study by Giampaolo and Illuminati Q yields instead that 
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there exist specific single-party unitary operations (corre- 
sponding to the projection on the z component of the spin 
for qubits and qutrits) which have the following proper- 
ties: (i) they leave a pure bipartite state invariant if and 
only if it is a product state; (ii) they transform any pure 
bipartite entangled state in such a way that the min- 
imum squared Euclidean distance of the original state 
from the set of all possible transformed states is an en- 
tanglement monotone. In the case of bipartite states of 
qubits and qutrits, such a measure coincides exactly with 
the marginal linear entropy, quantifyingthe degree of im- 
purity (mixedness) of each subsystem [8(. Therefore, en- 
tanglement monotones based on completely different def- 
initions, such as the linear entropy and the tangle [l(| , are 
re-discovered and re-interpreted in terms of the Hilbert- 
space distance between quantum states and their images 
under suitably selected local unitary operations. 

In this work, we apply the framework introduced in 
Ref. [1] to characterize entanglement of pure Gaussian 
states of continuous-variable systems. Recent progresses 
have showed that many nontrivial problems in entangle- 
ment theory, whose remarkable complexity renders their 
solution unachievable in quciit systems with d greater 
than 2 or 3, can be successfully tackled with different 
techniques when considering systems defined on infinite- 
dimensional Hilbert spaces, like, e.g. the quantum elec- 
tromagnetic field [Til] - In particular, Gaussian states, 
such as coherent, squeezed states, and in general all 
ground and thermal states of harmonic lattices, have 
played an increasingly important role in quantum infor- 
mation science, thanks to their simple structural proper- 
ties as well as to the high degree of experimental control 
on their production and manipulation [l2l . [T3 |. Moti- 
vated by these considerations, we seek here to provide 
a novel geometric interpretation for bipartite entangle- 
ment of pure Gaussian states, in terms of the perturba- 
tion induced on them by single-mode unitary operations 
in Hilbert space, or, equivalently, symplectic transforma- 
tions in quantum phase space. We will find, in direct 
analogy with the discrete- variable case 0] , that there ex- 
ists a single-mode symplectic operation which preserves 
product states, and whose action leads in general to 
the definition of a pure-state entanglement monotone for 
1 x N Gaussian states. This measure does not exactly 
coincide with any known entanglement measure, even 
though it is a monotonically increasing function of the 
entropy of entanglement, providing thus a novel quanti- 
fier of continuous variable entanglement endowed with a 
purely geometric interpretation. 

The paper is organized as follows: in Section II we 
briefly review the basic tools of the symplectic formal- 
ism in phase space, that is best suited for the analysis of 
separability, entanglement, and quantum operations on 
Gaussian states of infinite-dimensional quantum systems. 
In section III we introduce and analyze the properties 
of single-mode unitary (symplectic) operations in quan- 
tum phase space, and define the distance, induced by the 
fidelity, between pure Gaussian states and their images 



under such operations (these images are again pure Gaus- 
sian states). We then proceed to determine the minimum 
distance over the set of all possible such transformations, 
and the associated extremal operation. We prove that 
invariance of a state under the action of the extremal op- 
eration is a necessary and sufficient condition for the full 
separability of multimodc pure Gaussian states of trans- 
lationally invariant system, and show that the associated 
minimum distance is an entanglement monotone closely 
related to the linear entropy of the subsystem reductions. 
We finally discuss the relation between this novel entan- 
glement monotone and the various possible extensions of 
the definition of the tangle to continuous variable sys- 
tems. In Section IV we point out at some possible future 
lines of investigation in the framework of the formalism of 
local symplectic operations, also concerning mixed states, 
and discuss possible methods for the direct experimental 
detection of the minimum distance using linear optical 
elements, with an explicit example focused on tripartite 
Gaussian states. 



II. PHASE-SPACE DESCRIPTION OF 
GAUSSIAN STATES AND SINGLE-MODE 
SYMPLECTIC OPERATIONS 

We consider a continuous-variable (CV) system con- 
sisting of N canonical bosonic modes, associated with an 
infinite-dimensional Hilbert space, tensor product of the 
N single-mode Fock spaces [TJ [l2|, E3] • Unitary opera- 
tions which are at most quadratic in the canonical op- 
erators, amount to symplectic transformations in phase 
space. A real 2N x 2N matrix describes a symplectic 
transformation S € Sp(2N, R) if, by definition, it pre- 
serves the symplectic form, 

sns T = n, o = w ffiJV , u = ( _° x J ) . (i) 

For a single mode, the generators of the symplectic group 
Sp(2,R) are [H 




(2) 

where £2 = The matrices S,'s in Eq. @ arc traceless. 
Together with the identity matrix 1, they form a basis in 
the space of 2 x 2 real matrices. According to the Euler 
decomposition, the most general single-mode symplectic 
operation S € Sp(2, R) can be written as a sequence 
of a rotation, a squeezing, and a second rotation (with 
different angle) in phase space, 

g _ ( coscj) sin$> \ / £ \ / cos6> sin 9 \ , , 
— I — sin (j> cos (j> J I J \ — sin 9 cos 9 J ' ^ ' 

reducing to the identity transformation for 6 = <f> = 

(U = i. _ 

We are interested in studying the minimal distance be- 
tween a state and its image as transformed by a specific 
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type of local single-mode symplcctic operations. Clearly, 
one cannot allow the identity transformation in the defin- 
ing set of possible operations, if one wants to avoid end- 
ing up with a trivial null distance on all quantum states. 
Then, in analogy with the finite-dimensional case, we im- 
pose the condition of tracelessness Q, and we define a 
unitary single-mode operation S smo as the most general 
Sp(2, R) symplcctic transformation of the form Eq. 
with Tr S = 0. In this way we are only considering sym- 
plectic transformations which are orthogonal to the iden- 
tity. Imposing such constraint yields <f> = tt/2 — 9, namely 



(« 2 "i) 



cos x) sin f 



Ssmo — 



cos" sin 

1 



coA+sin 2 9 (g 2 -lfcosflsine I ' ( 4 ) 



The transformation S smo (£, 8) can be written as a linear 
combination of the E,'s from Eq. @, S smo = aT,i+f3T l2 + 
7E3, where the symplectic condition Eq. ([T]) imposes 7 = 
^/3 2 - a 2 - 1, (3 > Va 2 + 1. Explicitly [Hj: 



A pure Gaussian state is characterized by = 1, 
V k = 1...N, which implies Deter = 1. Such a state 
may be, for instance, the ground state of some harmonic 
Hamiltonian. We want to study the 1 x (TV — 1) entangle- 
ment of one mode with the remaining TV — 1 modes, via 
the perturbation induced by single-mode operations on 
mode 1. Namely, we aim to study the minimal squared 
distance between the Gaussian state <x and the state ob- 
tained from it by applying a S smo of the form Eq. ([5]) 
on any selected mode, for instance mode 1. It is impor- 
tant to recall that the transformed state, being obtained 
from the original pure Gaussian state by applying to it a 
symplectic transformation, i.e. a unitary transformation 
at most quadratic in the field variables, is again a pure 
Gaussian state. Introducing the standard Bures metric, 
the minimum distance reads 



D(a) 



min [1 

a,0 



T(cr, 



(7) 



Ssmo{d,(3) 



V/3 2 -a 2 -I 
a — (3 



a + (3 



-VW 2 



1 



(5) 



where the parameters a, (3 are connected with the squeez- 
ing £ and the rotation angle 8, appearing in Eq. ([5]), by 
the following relations: 



Here a' = [S smo (a, (3) © 1 2 ...n] ■ <* ■ [S smo (a, (3) © 1 2 ...jv] t , 
and the fidelity T between two pure-state TV-mode Gaus- 
sian CMs can be computed as |17l | 



T{a,a') = 2 N /^/Det {tr + tr'), 



i = \fia + ay/p-l]/c 



cos 9 



V>-oV^I-l]/[2(/3 2 -l)] 



(G) 



III. EXTREMAL SINGLE-MODE OPERATIONS 
AND ENTANGLEMENT OF PURE GAUSSIAN 
STATES 

We can now move to the specific setting of the ge- 
ometric analysis. Let our TV-mode bosonic system be 
prepared in a pure Gaussian state [ll|. We recall 
that Gaussian states of TV modes are completely de- 
scribed in phase space (once the first moments are set 
to zero via local displacements) by the real, symmet- 
ric 2TV x 2TV covariance matrix (CM) <r, whose entries 
are a lm = l/2({X u X m }) - (X ; )(l m )- Here X = 
{xi,pi, . . . , xn,Pn} is the vector of the field quadra- 
ture operators, whose canonical commutation relations 
can be expressed in matrix form: [X;,TV rn ] = 2if2; rn , 
with the symplectic form defined in Eq. |T]). Accord- 
ing to Williamson theorem [l6|, the CM of a TV-mode 
Gaussian state can be always diagonalizcd by means of 
a global symplcctic transformation (this corresponds to 
the normal mode decomposition): W^crW^ = where 
Wo- £ 5p(2TV,R) and u = 0f =1 diag{i/ fe , v k } is the CM 
corresponding to the tensor product of single-mode ther- 
mal states. The quantities {vk] are the so-called sym- 
plectic eigenvalues of the CM cr. 



amounting to the overlap |('0|'0 / )| 2 between the original 
and the perturbed Gaussian quantum states. 

To proceed in the evaluation of Eq. 0, let us first 
take into account that, in full generality, pure Gaussian 
states can always be brought in the phase-space Schmidt 
form 11] with respect to the 1 x (TV — 1) bipartition. 
The symplectic transformation W achieving the Schmidt 
decomposition is the direct sum of the two Williamson 
diagonalizing operations acting on the single-mode and 
the (TV — l)-mode subspaces, respectively, W = W\ © 
W2...N • One then has 



cr w 



WerW 1 



a 


-VaJ~7J 




a 




v a 







a 



L3...JV , 



with a > 1 , 



(8) 



i.e. the phase-space Schmidt form of a is constituted by 
one two-mode squeezed state between modes 1 and 2, 
tensor TV— 2 uncorrelated vacua [l8[. To evaluate Eq. (J7J), 
we need the minimum of Det [<T+(S S mo®^-2...N)c(S sm o® 
^2...n)] T ■ We will now show that it is enough to consider 
states in the form cry/- In fact, 
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Det [<rw + {Ssmo © ^-2...n)<^w{S srno © 3-2.. jv) ] 
= Dot [(Wi © W 2 ...n)(t(W? © Wl jf) + {S smo Wi © W 2 ...n)(t(W? Sj mo © W 2 T _ N )] 
= Det [(Wi © W 2 ... n )ct(W[ © Wl, N ) + {WxW^SsmoWx © W 2 ... n )<t(W[ Sj mo W^ 1T W[ © W 2 t .jv)] 
= Det{[Wi © W 2 ...jv][o- + (Wf ^oWi^.jvMWf l S am0 Wi © 3 2 ... w ) T ][Wi © IF 2 ..jv] t } 
= Dot [a + {W- 1 S smo W 1 © l 2 ... N )(r(Wr 1 S smo W 1 © 1 2 ...tv) T ] , (9) 



where we exploited the group properties of Sp(2, R), the 
fact that a symplectic operation S has Det 5 = 1, and 
the property that the inverse S^ 1 of a symplectic trans- 
formation S is itself symplectic. Now, from the cyclic 
property of the trace, it follows that W^ 1 S smo Wi is it- 
self a traceless symplectic operation, i.e. a single-mode 
operation of the form Eq. (J5j) . Thus the minimum of the 
above determinant, taken over the entire set S smo (a, f3) 
of single-mode unitary operations, is invariant under lo- 
cal symplectic operations Wi© W 2 ...n performed on state 
<r. Thus, without loss of generality, we can choose a pure 
iV-mode Gaussian state in the phase-space Schmidt form 
cr w of Eq. (0. Therefore, 



Det [ar w + {Ssmo © 1 2 ...n)<tw(S s 

= 2 2 < Ar - 1 )[( a 2 -l) 2 + 4a 2 /3 2 ]. 



12...JVJ 



(10) 



The minimum is then acquired, as fj > va 2 + l, for 
/3 = 1, a = 0. The corresponding extremal single- mode 
operation is then, finally 



S sm o{0 : 1) 



1 

-1 



(11) 



This is a simple rotation of tt/2 in phase space, and may 
be seen as the CV analogue of the spin-flip operation on 
qubits, realized by the a y Pauli matrix. 

We observe that a product state, characterized by a 
CM in direct sum form, er® = o"i©<72..jv, is left invariant 
by the extremal single-mode operation: 



[S smo (0, 1) © l 2 ...Ar]o- ffi [S smo {0, 1) © 1 2 . 



.N 



(12) 



That is, on pure product Gaussian states, extremal and 
invariant (or preserving) operations coincide, in full anal- 
ogy with the finite-dimensional case analyzed in Ref. Q . 
Hence, a pure Gaussian states is separable if and only if 
there exists a traceless single-mode symplectic (unitary) 
operation that leaves it unperturbed. This is again in per- 
fect analogy with the discrete- variable analysis performed 
for qubits and qutrits 0]. The minimum distance Eq. |7]). 
achieved for /3 = 1, a = 0, can now be evaluated explic- 
itly and reads 



D(a) = 1 



1 



2 JV " 1 (a 2 + l) a 2 + l 



(13) 



The quantity D(cr) is a measure of the entanglement 
between mode 1 and the rest of the system, being an 



increasing function of the single-mode mixedness factor 
a > 1. For product states a = 1 and one correctly re- 
trieve D{cr®) = 0. One should recall that Eq. (H3]) holds 
in general, and not only for states in Schmidt form, once a 
is identified with the (locally invariant) symplectic eigen- 
value of the reduced CM <j\ of mode 1, computable from 
its determinant as a — VDct cr\. From this observation, 
it follows immediately that D(cr), even if constructed 
through the action of local unitaries, is invariant under 
them - as already proved by Eq. (0) - and is thus a proper 
entanglement measure. 

The 1 x (N — 1) linear entropy of entanglement for the 
state er (corresponding to the tangle for qubits [lj|) reads 



E l (<t) 



1 



(14) 



We see that El is a monotonic function of D, thus 
qualitatively equivalent to but yet not exactly coincident 
with the latter everywhere, at variance with the discrete- 
variable case, in which they do strictly coincide [|[ (The 
behavior of the two entanglement monotones is compared 
in Fig. [1]). The fact that the linear entropy of the re- 
duced state does not coincide exactly with the minimum 
distance achieved under local symplectic operations may 
be traced back to the non uniqueness in the definition 
of the "tangle" for Gaussian states of CV systems. For 
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FIG. 1: (color online) Entanglement of pure (1 X (N — 1))- 
mode Gaussian states as a function of the single-mode squeez- 
ing r, defined such that a = cosh 2r. The solid line depicts 
the distance-based geometric measure D defined in Eq. (|13[1 . 
while the dashed line corresponds to the linear entropy of 
entanglement El, Eq. {H}. 
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qubits, at least four different definitions coalesce into the 
same entanglement monotone: (i) squared concurrence 
f3~ot | ; (ii) local linear entropy [l9| ; (iii) squared negativity 
(negativity equals concurrence for pure qubit states [20l|): 
(iv) minimum distance under single-qubit unitary trans- 
formation Q . On the other hand, while the concurrence 
is not well defined in CV systems, the other definitions of 
the tangle all give rise to different (yet equivalent) entan- 
glement quantifiers in these systems. For instance, the 
Gaussian tangle defined as the squared negativity [21j , in 
analogy with definition (iii), reads 



T G (cr) = [a 2 + aV 'a 2 - 1 - l]/2 . 



(15) 



The von Neumann entropy of entanglement, for refer- 
ence, is given by E v (ct) = [{a + l)/2]log[(a + l)/2] - 
[(a — l)/2] log[(a — l)/2]. All these measures are mono- 
tonically increasing functions of each other (and of a), 
some of them being normalized between and 1 (like 
D and -Ex), the others diverging in the limit of infinite 
squeezing, a — > oo. 



IV. EXPERIMENTAL REMARKS AND 
FUTURE PERSPECTIVES 

The minimum distance D provides a new bona fide geo- 
metric measure of entanglement for pure Gaussian states, 
close in spirit to the low-dimensional, discrete-variable 
counterpart introduced in Ref. ,8|, and embodying yet 
another generalization of the tangle. However, we would 
like to remark that, among the three possible CV versions 
of the tangle, only tq, Eq. (fl5|) . satisfies the CV gener- 
alization of the Coffman-Kundu-Wootters monogamy in- 
equality [13, E^Ij as proved in Ref. [2l[ for all, pure and 
mixed, V-mode Gaussian states. 

On the other hand, the geometric measure of entangle- 
ment D(cr) that we have introduced in this work for pure 
Gaussian states, has the nice property of being amenable 
to direct experimental investigation, once two copies of 
an unknown Gaussian state with CM er are available. 
One first needs a phase plate in order to rotate one copy 
of 7r/2, realizing the operation £2, as demonstrated e.g. in 
f22| . Thereafter, the evaluation of the overlap between 
the rotated copy and the unrotated one involves stan- 
dard tools of linear optics, as routinely demonstrated in 
the determination of the fidelity T of teleportation ex- 
periments with continuous variables [23| J or in the im- 
plementation of intcrfcromctric schemes [24| that can be 
realized even without homodyning [25| . Our result thus 
suggests a way to the direct measurement of CV entan- 
glement in pure 1 x N Gaussian states, in analogy with 
what achieved experimentally in the case of qubits: In 
that case, the entanglement, quantified by the two-point 
concurrence, has been directly measured on the two-fold 
copy of unknown two-qubit pure states [26| . 

This proposal looks especially appealing for Gaussian 
states with a small number of modes. A relevant example 
is provided by three-mode Gaussian states, whose CM 




assumes in general the following expression in terms of 2 
by 2 submatrices, 



(16) 



The structural and informational properties of three- 
mode Gaussian states, with a special emphasis on the 
pure-state instance, have been extensively studied in 
Ref. p7j . while a scheme for their production via in- 
terlinked nonlinear interactions in x' 2 ' media has been 
presented in Ref. [28[, together with preliminary exper- 
imental results. When modes 2 and 3 have the same 
average number n of thermal photons, the correspond- 
ing (parametric) pure three-mode Gaussian state is said 
to be 'bisymmetric' and its CM can be written in the 
standard form of Eq. (fT6|) , with 



0"2 = CT3 

£23 

£12 = £13 



al2, 
a+ 1 
2 

a- 1 



1 2 , 
1 2 , 



diag ■ 



Va 2 ^! Va 2 ^! 



V2 



V2 



and a = An + 1. 



The geometric entanglement between the first mode and 
the group of modes 2 and 3, as obtainable from the single- 
mode perturbation Eq. (fTTj) applied on mode 1, is then 
directly given by Eq. (fl"3"|) as a function of n. The three- 
mode Gaussian states of this family are known to be op- 
timal resources for 1 — > 2 CV telecloning (i.e. cloning at 
distance, or equivalently teleportation to more than one 
receiver) of sing le-mode coherent states [29| , as discussed 
also in [2?], l28j . The single-clone fidelity T exhibits a 
non-monotonic, concave behaviour as a function of h, 
reaching the maximum jF max = 2/3 at n = 1/2. Very 
recently, the first experimental demonstration of uncon- 
ditional 1 — > 2 telecloning of unknown coherent states has 
been realized [3(| , with a measured fidelity for each clone 
of T = 0.58 ± 0.01 (surpassing the classical threshold of 
0.5). This experimental milestone has raised renewed in- 
terest towards CV quantum communication [3l| . In the 
context of this work, such an achievement entails that all 
the elementary steps required to access pure-state Gaus- 
sian entanglement from a geometric point of view have 
been already successfully undertaken. Our prescription, 
therefore, is likely to be seen "at work" experimentally 
on multimode Gaussian states in the near future. 

In this paper we have dealt with pure Gaussian states 
only. It is natural to ask whether a suitable generaliza- 
tion of the present approach is able to provide a geometric 
interpretation, possibly amenable to direct experimental 
tests, of mixed-state entanglement measures as well. In 
this respect, it is important to clarify to which extent the 
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results of this paper are still valid for mixed states. In 
primis, it is generally true that the extremal single-mode 
operation, Eq. (fTTj). preserves tensor product, even mixed 
Gaussian states [see Eq. (|12p]. However, convex combi- 
nation of product states, i.e. separable mixed Gaussian 
states, are not left invariant by the action of such local op- 
eration. Accordingly, the corresponding geometric mea- 
sure (minimum distance D) defined by Eq. (|13|) is not, in 
general, an entanglement monotone for mixed Gaussian 
states. One can thus conclude that, in the mixed-state 
scenario, the mere action of S smo (Q, 1) leads to a distinc- 
tion between tensor product states (totally uncorrelated, 
on which the distance Eq. (JT3J) is zero) and states display- 
ing some form of (quantum and/or classical) correlation. 
A refinement is henceforth necessary in order to aim at 
discriminating, from a geometric point of view, the quan- 
tum portion entanglement - from the total amount of 
correlations. A feasible way to deal with this issue seems 
that of identifying a minimal set of single-mode unitary 
operations, such that a suitably defined "distance" in- 
volving their combined action, may turn to be equivalent 
(or to provide bounds) to known entanglement mono- 
tones (e.g. negativities, tangles and/or measures based 



on the Gaussian convex roof [32(1). One should then be 
able to readily provide a recipe for the practical estima- 
tion of mixed-state Gaussian entanglement with few local 
measurements (see also [33^. 

Finally, we would like to remark that the framework in- 
troduced in Ref. Q , and further discussed in the present 
paper, can be naturally applied to investigate criticality 
and entropy scaling in the ground states of harmonic lat- 
tices [Hj], with the purpose of establishing connections 
similar to those unveiled for the ground states of spin 
systems at criticality |3a |. This interesting perspective 
will be the object of further future studies. 
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